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ABSTRACT 

We have discovered a detached pair of white dwarfs (WDs) with a 12.75 min orbital period and a 
1,315 km s~^ radial velocity amplitude. We measure the full orbital parameters of the system using its 
light curve, which shows ellipsoidal variations, Doppler boosting, and primary and secondary eclipses. 
The primary is a 0.25 M0 tidally distorted helium WD, only the second tidally distorted WD known. 
The unseen secondary is a 0.55 Mq carbon-oxygen WD. The two WDs will come into contact in 0.9 
Myr due to loss of energy and angular momentum via gravitational wave radiation. Upon contact 
the systems may merge yielding a rapidly spinning massive WD, form a stable interacting binary, or 
possibly explode as an underluminous supernova type la. The system currently has a gravitational 
wave strain of 10~^^, about 10,000 times larger than the Hulse- Taylor pulsar; this system would be 
detected by the proposed LISA gravitational wave mission in the first week of operation. This system's 
rapid change in orbital period will provide a fundamental test of general relativity. 
Subject headings: binaries: close — binaries: eclipsing — Gravitational waves — Stars: individual: 
SDSS J065133. 338-1-284423. 37 — white dwarfs 



1. INTRODUCTION 

In Einstein's general theory of relativity, close pairs of 
stars produce gravi tational waves , ripple s in the curva- 
ture of space-time (jEinsteinl 119161 11918[) . Observations 
of the Hulse- Taylor binary pulsar PSR B1913+16 con- 
firm the slow decay of the orbit predicted by Einstein's 
theory and provide iiidirect evidence for g ravitational 
waves pulse fc Tavloilll975tlWeisberg et al.l[2010) . On- 
going and proposed ins truments seek to detect gravi- 
tational waves directly (pafrv et aTl 119941 : iHobbs et al.l 
I2009H Abbott et al.|[2010t ). The strongest known gravita- 
tional wave sources are compact binaries containing neu- 
tron stars and WDs (e.g., |Nclemans 2009). There are 
presently four known binaries with orb ital periods less 
than 15 minutes - three AM CVn stars (Habe rl fc MotchI 
ri995l: llsrael eFaniTggl IWarner fc Woudt 200D and one 
low mass X-ray binary (jStella et al.iil987i) - but all bi- 
naries with < 1 hr orbital periods, with one exception 
(jKilic et al.ll2011 m , are interacting s ystems that compli- 
cate d irect tests of Einstein's theory. iMarsh fc Nelemanj 
(|2005[ ) show that the expected period change due 
to gravitational wave radiation can be modified by 
both induction-driven angular momentum interchange 
in magnetic WD systems and by mass-transfer in ac- 
creting systems. Moreover, there are different ex- 
planations for both the nature of the shortest-period 
AM CVn systems an d their observed period changes 
JStrohmavcr 2004; D'A ntona et al.ll2006HDelove fc Taaiiil 
[2006t iRoelofs et al.ii2010[l 
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Here, we report the discovery of a detached, 765 sec 
orbital period binary, SDSS J065133.33+284423.3 (here- 
after J0651). The discovery comes from our on-going 
ELM Survey, a targeted spectrosc opic survey for ex- 
tremely low ma ss <0.25Mq WDs (jBrown et al.l 120101 : 
iKilic et al.ll2010D . The ELM Survey demonstrates that 
extremely low mass WDs are found in compact bi- 
naries, and that the majority of these binaries have 
<10 Gyr merger times due to gravitational wave radi- 
ation (Kilic et al. 2Glla). Depending on the poorly con- 
strained physics of the merger process, these merging 
WD systems may produce single helium-rich sdO stars, 
stable mass-transfer AM CVn binarie s, or possibly un - 
derluminous supernova (as discussed in lKilic et al.ll2010() . 
The observed merger rate of the low mass WD binaries is 
about 1% of the Type la supernovae rate - comparable 
to th e rate of underluminous supernovae (|Brown et al.l 
l2011f ). The ELM Survey is a lso responsible for finding 
the first tidally distorted WD (jKilic et al.ll2011bD . J0651 
is the most extreme system discovered to date and shows 
that the ELM Survey may open a new window on grav- 
itational wave sources. 

In J0651, the orientation of the binary allows us to 
observe eclipses of each star by the other, leading to ac- 
curate measurement of the orbital parameters, masses, 
and WD radii. There is no evidence for mass transfer. 
These results suggest that J0651 is the cleanest known 
strong gravitational wave source, with an orbital decay 
time of less than 1 Myr. Future observations will allow us 
to measure the change in the orbital period. We hope to 
one day compare this change with direct measurements 
of gravitational waves and provide an unprecedented test 
of general relativity. 
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Fig. 1. — The composite observed spectrum of J0651, con- 
structed by shifting the individual spectra to rest-frame and then 
summing them together to create an effective 60.5 min exposure, 
compared to the best-fit stellar atmosphere (solid red line) for a 
0.25 Mq WD. The broad hydrogen Balmer absorption lines visu- 
ally indicate J0651 is an unusually low mass WD. 



2. DATA AND ANALYSIS 

We discovered J0651 on 2011 March 2 as part of our 
targeted spectroscopic survey for extremely low mass 
WDs using the 6.5m MMT telescope. We used the MMT 
Blue Channel spectrograph (Schmidt et al. 1989) to ob- 
tain 1 A resolution spectra in the wavelength range 3600 
A - 4500 A. We recognized J0651 as a low mass WD due 
to its pressure-broadened hydrogen Balmer lines, seen in 
Figure [TJ 

Spectral fits to hydrogen atmosphere WD models 
(|Koestei1l2008D yield a surface gravity of log 5 = 6.79 ± 
0.04 dex {g in cm s~^) and an effective temperature of 
Teff = 16, 400 ± 300 K. T he Sloan Digital Sky Survey 
(SDSS, lAihara et al.ll20ilT ) spectrum of this object, on 
the other hand, yields a systematically larger logy = 
6.97 ± 0.05 dex and Teff = 17, 700 ± 250 K (S. J. Klein- 
man 2011, private communication) because the SDSS ex- 
posure spans 3.5 orbital periods and thus measures arti- 
ficially broadened Balmer lines. Our own 6-10 min ex- 
posures of J0651 show similarly over-estimated surface 
gravities and temperatures. A correct understanding of 
J0651 requires time- resolved spectroscopy only possible 
with a larger telescope like the 6.5m MMT. 

Figure [2] compares our best-fit model atmosphere to 
dereddened {E{B-V) = 0.0706 mag) ultraviolet Galaxy 
Evolu tion Explorer (GALEX) photometry (|Martin et al.l 
l2005f ). optical SDSS photometry, and our own near- 
infrared J = 1 9.599 lb 0.029 mag m easurement obtained 
with SWIRC ()Brown et al.l [2008[) on the MMT. The 
broadband photometry supports our spectroscopic effec- 
tive temperature and surface gravity measurement with 
one possible exception. The GALEX NUV point shows 
a 2-CT discrepancy with the model spectrum convolved 
with the filter bandpass, however the photometric error 
does not include the significant uncer tainty in the extinc- 
tion correction for GALEX fiUers (jWvder et al.l I2005D 
and so we consider the discrepancy suspect. Compari- 
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Fig. 2. — Best-fit WD model atmosphere (dotted line) compared 
to broadband photometry (dots). The ultraviolet, optical, and 
near-infrared measurements support our spectroscopic fit. 



son with WD evolutionary tracks (jPanei et al.l 120071 ) in- 
dicates that J0651 is consistent with an extremely low 
mass 0.25 M0 helium-core WD. 

We discovered that J0651 is a compact binary system 
when back-to-back spectra separated by 6 min showed 
a ~1300 km s~^ change in radial velocity (Figure [3]). 
Despite the faintness of the g = 19.1 WD, the light- 
gathering power of the 6.5m MMT telescope allowed us 
to reduce our exposure times to 2-2.5 minutes and re- 
solve the orbital period. The data are presented in Table 
[TJ The observed velocity amplitude is an underestimate, 
however, because 1) our exposures span 18% of orbital 
phase and 2) the radial velocity curve is sinusoidal, not 
linear. By integrating a sine curve at the phase covered 
by our exposures, we determine that the velocity ampli- 
tude correction is 5.5%. The observed velocity amplitude 
is also underestimated if the WD is tidally distorted and 
its center-of-light differs from its center-of-mass. We do 
not correct for this effect given that the velocity correc- 
tion is less than 1% - comparable to our measurement 
uncertainty - for the observed oblateness. 

The best-fit orbital parameters to the radial veloci- 
ties are: period P = 765.4 ± 7.9 sec, corrected velocity 
semi-amplitude K = 657.3 ± 2.4 km s~^, and systemic 
velocity 7 = 16.6 ± 0.6 km s~^. J0651 is very likely 
a Galactic disk object based on i ts small proper mo tion 
{fia, fJ-s) = (-3.6, -1.2) mas yr"! (jMunn et al.ll2004D and 
small systemic velocity. J0651's location 220 pc above 
the plane (for a distance of 1 kpc, see below) is also con- 
sistent with a disk object. 

We obtained time-series optical photometry of J0651 
using the McDonald 2.1m Otto St ruve Telesc ope using 
the A rgos frame transfer camera (jNather fc: Mukadaml 
|2004( ). Images were obtained with a BG4O Schott glass 
filter. Our set of 3,038 10- and 15-second exposures come 
from five different nights in 2011 April covering a total 
time baseline of 12.1 days. The observed light curve, 
phased to the best-fit period, is plotted in Figure S] and 
shows three significant features: a sinusoidal pattern due 
to ellipsoidal variations from the tidally distorted WD, 
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Fig. 3. — Radial velocity observations phased to the 765 sec 
orbital period. The best-fit orbit (dotted line) has a 1314.6 km 
s~^ velocity amplitude, wrhich is 0.44% the speed of light. 



TABLE 1 
Radial Velocity Measurements 



HJD+2455600 


Vhelio 


(days)t 


(km s-1) 


22.76724617 


471.0 ± 30.2 


22.76919045 


-310.7 ± 28.9 


22.77118102 


-534.9 ± 27.8 


22.77311373 


189.3 ± 23.9 


22.77506958 


684.5 ± 23.1 


22.77697914 


132.2 ± 25.0 


22.77892341 


-555.3 ± 38.4 


22.78221017 


177.1 ± 38.3 


22.78415444 


568.7 ± 29.0 


22.78607558 


34.5 ± 27.7 


22.78806615 


-505.9 ± 41.3 


22.79001043 


-172.5 ± 25.4 


22.79194313 


540.2 ± 28.0 


23.66752771 


92.8 ± 24.4 


23.66963400 


537.0 ± 22.3 


23.67142783 


140.8 ± 29.7 


23.67304806 


-544.2 ± 34.5 


23.67469144 


-593.0 ± 26.1 


23.67628852 


218.9 ± 44.2 


23.67788561 


496.0 ± 29.6 


23.68036224 


124.7 ± 25.4 


23.68195933 


-466.8 ± 25.0 


23.68353327 


-580.2 ± 27.4 


23.68510721 


26.8 ± 22.8 


23.68670429 


669.3 ± 25.5 


23.68830137 


471.9 ± 26.0 


23.68992161 


-193.7 ± 24.5 



Note. — fBased on UTC. 

an asymmetric peak in light due to relativistic beaming 
(so-called Doppler boosting), and periodic dips in light 
from the eclipses of the primary and secondary WDs. 

We model the light curve of J0651 using JKTE- 
BOP (|Southworth et all [200l) and verify our resuhs 
with PHOEBE (jPrsa fc ZwitteH[2005l) . JKTEBOPQand 

^ JKTEBOP models the projection of each star as a biaxial 
ellipsoid and calculates a light curve by numerical integration of 
concentric circles over each star. A by-product of this calcula- 



FlG. 4. — J0651 light curve. The upper panel plots the observed 
photometry vs. orbital phase, while the lower panel compares the 
binned data to the best-fit model (solid red line). The data reveals 
three dramatic features: a sinusoidal pattern due to ellipsoidal 
variations from the tidally distorted WD, an asymmetric peak in 
light due to Doppler boosting, and periodic dips in light from the 
eclipses of the primary (at phase 0) and secondary (at phase 0.5) 
WDs. 

PHO EBE are based on the Eclipsin g Binary Orbit Pro- 
gram dPopper fc Etzel|[l981f ) and the lWilson fc DevinnevI 
(1971) codes, respectively. We first remove the 0.5% 
Doppler boosting (relativistic beaming) signal, how- 
ever, because neither code models Doppler boosting. 
Doppl er boosting has been seen i n only a handful of sys- 
tems (Ivan Kerk wiik et"aD ' 20101: iMazeh fc Faigle^ 120101: 
iShporer et al.l l2010: Blocme n et al.ll201lD : its asvmmet- 
ric contribution to the J0651 light curve confirms the 765 
sec orbital period. We then fit the 5% amplitude ellip- 
soidal variations and the 15% primary eclipse and 4% 
secondary eclipse. Reflection effects due to the heating 
of each WD by its companion are weak, 0.3%±0.4%. 

The light curve yields a much more precise measure- 
ment of orbital period, P = 765.2062 ± 0.003 sec. The 
ellipsoidal variations are due to the changing projected 
area of the distorted primary WD, hence they strongly 
depend on the inclination angle of the binary system. 
J0651's ellipsoidal variations and eclipses constrairlJ the 
orbital inclination to i = 86.91;]^q deg. 

Eclipses also provide a precise measurement of the WD 
radii. The 0.25 Mq primary WD has an observed radius 
of 0.0353±0.0004 R© that differs by 5% fro m the 0.0337 
Rf:) r adius predicted by helium WD models (jPanei et alj 
120071 ). Going in the other direction, the models predict 
that a WD with the observed radius and with mass 0.24 
Mq has log 5 = 6.71 dex and Toff = 16, 000 K, consistent 
with our spectroscopic observations. Given the uncer- 
tainties, we consider the observed and predicted radii in 
excellent agreement. 

The unseen secondary has a mass of 0.55 Mq and 
an observed radius of 0.0132±0.0003 R© typical for a 

tion is an estimate of the oblateness of the primary WD. The 3% 
oblateness of J0651 falls within JKTEBOP's 4% oblateness limit 
for accurate light curve analysis. 

* Errors are estimated fr om 10,000 Monte-Carlo si mulations with 
JKTEBOP as described in lSouthworth et al.l ifSOOSD . 
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carbon-oxygen WD. The eclipse depths indicate the sec- 
ondary is 3.7% ± 0.3% as bright as the primary. Adopt- 
ing Mg = 8.9 ± 0.1 mag a s the a,bsolut e magnitude of 
the 0.25 M0 primary (Pane iet al.ll2007l ). the secondary 
has Mg = 12.5 mag. A carbon-oxygen WD of this lu- 
minosity has Toff = 9000 K, \ogg = 7.9 dex, a cooling 
age of 700 Myr, and a radius of 0.0137 R© (P.-E. Trem- 
blay 2011, private communication). The observed sec- 
ondary radius differs by 4% from the predicted radius. 
There arc only two other model-independent mass and 
radius determinations of helium and carbon -oxygen WDs 
(Steinfadt et al.ll2010t iParsons et al]l2011[ ). Our results 
thus provide important constraints and reveal an overall 
agreement with current WD models. 

3. DISCUSSION 

Having accurately measured J0651's binary param- 
eters we now calculate the gravitational wave emis- 
sion predicted by Einstein's general theory of relativ- 
ity. J0651, with a de-reddened apparent magnitude 
go = 18.84 ± 0.01 and thus a 1.0 ± 0.1 kpc distance 
from the Sun, has a predicted gravitational wave strain 
(e.g.. IR,oelofs"ei~aLi 20071 oi h = 1.2 x 10-^2. J0651 is 
among the strongest known gravitational wave sources 
and, more importantly, has an orbital frequency that 
places it well above the ex pected gravitational wave fore- 
ground (Nelemans"2009'). The proposed ES A/NASA 
Laser Interferometer Space Antenna (LISA) mission has 
its peak sensitivity at frequencies corresponding to ~5 
min orbital periods, and it should detect a strain of 10^^^ 
at th ese frequencies with a signal-to-noise of ~100 in one 
year (i Roelofs et al.|[2007l ). In other words, J0651 is a ver- 
ification source that LISA, if built, would detect in the 
first week of operation. 

We also predict that J0651's orbital period is shrink- 
ing by 2.7 X 10~^ sec per year due to gravitational wave 
radiation. The expected change in period adds up to a 
5.5 sec change in time-of-eclipse in one year. When we 
measure this change we expect to provide yet another 
fundamental test of general relativity and the existence 
of gravitational waves. 

The absence of mass transfer in J0651 is perhaps sur- 
prising given how quickly it will merge. In all known bi- 
naries with periods comparable to J0651 fe.g.. lNelemansl 
[2009) one star fills its Roche lobe and transfers mass to its 
companion. Our data show that the J0651 primary has a 
Roche lobe radius 1.5 times its current radius. Under the 
assumption of energy and angular momentum loss due to 
gravitational wave radiation, the primary WD will reach 
its Roche lobe radius at an orbital period of 6.5 minutes 
in 0.9 Myr. 

What happens when the primary WD fills its Roche 



lobe and begins mass transfer in 0.9 Myr is an open theo- 
retical question. The stability of the mass transfer phase 
depends in part on the known binary mass ratio q = 0.45 
and in part on the unknown tidal syn chronizati on and en- 
tropy of the primary WD (jMarsh et al.l 12004"). If mass 
transfer is stable, J0 651 will becorne an AM CVn sys- 
tem like HM Cancri (jRoelofs et al.ll2010[) . Recent WD 
mass transfer models, on the other hand, predict that 
mass transfer wil l cause the pair o f WDs to quickly coa- 
lesce and merge (jDan et al.ll2011l ). A merger will likely 
produce a rapidly spinning massive WD, but an underlu- 
minous supernova explosion is also a possibility. Under- 
luminous supernovae, such as SN 2005E, are rare types 
of supernovae that are 10-100 times less luminous than 
a normal super nova Type la and have only ~0.25 M© 
worth of ejecta (jPerets et al.ll2010l ). The spectrum and 
light curve of SN 2005E, for example, can be explained 
by the detona tion of a 0.2 Mg, helium layer on a 0.45 
Mq WD (Wal dman et al.ll2010[ ) - parameters very sim- 
ilar to the J0651 system. Depending on the nature of 
the mass transfer, there may be other mechanisms by 
which the system could detonate (jBildsten et al.l 120071 : 
Guillochon et al. 2010). Completing our survey for ex- 
tremely low mass WD systems such as J0651 will allow 
us to measure the space density and merger rate of these 
systems and compare with the rate of underluminous su- 
pernovae. 

4. CONCLUSION 

The eclipsing detached WD binary system J0651 
presents us with a remarkable laboratory. We can use 
its eclipses to measure WD masses and radii for detailed 
tests of WD models. We can use its changing binary 
orbital period and gravitational wave strain to test for 
the existence of gravitational waves predicted by general 
relativity. We can use its very existence to constrain the 
space density and merger rate of low mass WD binaries 
and links to underluminous supernovae. In the future, we 
plan to use multi-passband photometry to directly mea- 
sure the nature of the secondary WD and to detect the 
change in orbital period predicted by general relativity. 
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